First observation of PeV-energy neutrinos with IceCube 
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We report on the observation of two neutrino-induced events which have an estimated deposited 
energy in the IceCube detector of 1.04 ± 0.16 and 1.14 ± 0.17 PeV, respectively, the highest neu- 
trino energies observed so far. These events are consistent with fully contained particle showers 
induced by neutral-current v^.^.t (^e.n.r) or charged-current (i? e ) interactions within the IceCube 
detector. The events were discovered in a search for ultra-high energy neutrinos using data corre- 
sponding to 615.9 days effective livetime. The two neutrino events are observed over an expected 
atmospheric background of 0.082±0.004(stat)^g'g^(syst). The resulting p- value for the background- 
only hypothesis is 2.9 x 10 -3 (2.8<r) taking into account the uncertainty on the expected number of 
background events. Though the two events could be a first indication of an astrophysical neutrino 
flux, the moderate significance and the uncertainties on the expected atmospheric background from 
neutrinos produced in the decay of charmed mesons do not allow for a firm conclusion at this point. 



PACS numbers: 98.70.Sa, 95.55.Vj 

Astrophysical neutrinos are key probes of the high- 
energy universe. Because of their unique properties, 
neutrinos escape even dense regions, are undeflected in 
galactic or extra-galactic magnetic fields and traverse the 
photon- filled universe unhindered. Thus, neutrinos pro- 
vide direct information about the dynamics and interiors 
of extra-ordinary objects in deep space like gamma-ray 
bursts and active galactic nuclei. Neutrinos at energies 
above several 100 TeV are particularly interesting as the 
atmospheric background in this region is very low and a 
few astrophysical neutrinos can be significant. This let- 



ter reports on the observation of two high-energy particle 
shower events discovered in a search for ultra-high energy 
neutrinos above about 1 PeV using the IceCube detector. 

IceCube [l[ detects and reconstructs neutrinos by 
recording Chcrcnkov photons emitted from secondary 
charged particles produced in neutral-current (NC) or 
charged-current (CC) interactions of the neutrinos in the 
2,800 m thick glacial ice at the geographic South Pole. 
IceCube construction started in 2005 and was completed 
in December 2010. The detector consists of an array 
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FIG. 1. Surface view of the full IceCube (IC86) detector lay- 
out. Filled marks represent the x-y positions of the IceCube 
strings. Red marks in the central region are the DeepCore 
strings. Squares represent the strings that did not exist in 
the IC79 configuration. Open circles are the positions of the 
closest strings to the observed two cascade events. Stars are 
the reconstructed vertex positions of the two cascade events. 



of 5,160 optical sensors (digital optical modules, DOMs) 
on 86 strings at depths between 1,450 m and 2,450 m 
that instrument a volume of 1 km 3 of ice. Eight out of 
the 86 strings belong to the DeepCore sub-array [2j, a 
more densely instrumented volume in the bottom cen- 
ter of the detector. Each DOM consists of a 10" photo- 
multiplier tube (PMT) that is protected from the high 
pressure of the ice by a glass sphere that also contains 
electronics to amplify and digitize the pulses [3j from the 
PMTs |4J . Events are recorded as a series of pulses (wave- 
form) in each DOM where one can distinguish between 
two basic neutrino event signatures: a track-like light 
pattern in the detector originating from neutrino- induced 
muons (tracks) and a spherical light pattern produced by 
hadronic or electromagnetic particle showers (cascades). 
If not explicitly mentioned, we do not distinguish be- 
tween neutrinos and anti-neutrinos here. 

The analysis selects neutrino candidates calorimetri- 
cally using the total number of observed photo-electrons 
in each event (NPE) as a proxy of the deposited en- 
ergy thus retaining both bright tracks and cascades. 
Background events for this analysis are induced by muons 
and neutrinos generated in interactions of cosmic-rays in 
the atmosphere. Because of their steeply falling energy 
spectra, little background is expected in the signal re- 
gion above 1 PeV. The zenith angle distribution of at- 
mospheric muons peaks in the downward-going direction 
and sharply decreases towards the horizon with a cut-off 
at a zenith angle 8 of cos 9 ps 0.15 due to absorption in 
the Earth, while the atmospheric neutrino distributions 
have a weaker zenith-angle dependence. In addition to a 
global lower cut on NPE, the analysis rejects downward- 
going atmospheric muons by employing event reconstruc- 



tions based on a track hypothesis in combination with a 
higher NPE cut in the downward-going region. All re- 
maining events above the combined NPE threshold are 
considered to be signal candidates independent of their 
topological properties. 

Data were collected between May 2010 and May 2012, 
an effective livetime of 615.9 days excluding 54.2 days 
used for the optimization of the analysis. From May 
2010 to May 2011, DOMs on 79 strings (IC79) were oper- 
ational. This IC79 run period was immediately followed 
by the first year of data taking with the full 86-string 
(IC86) detector that lasted until May 2012. The layout 
of the IC86 configuration is shown in Fig. Q] Events are 
triggered when eight or more DOMs record signals in lo- 
cal coincidences. Such a coincidence occurs when a near- 
est or next-to-nearest DOM on the same string triggers 
within ±1 /is 

The collected data are filtered at the South Pole with a 
condition NPE > 1000, and then sent to a northern com- 
puter farm via satellite. In order to avoid biases, only 
about 10% of the data was used to develop the event se- 
lection. Photo-electron arrival times are extracted from 
each waveform and stored as "hits" . To remove random 
coincident noise hits, a two-staged hit cleaning based 
on the spatial distance and the time interval between 
hits is applied. Data from the DeepCore strings is dis- 
carded to maintain uniformity across the detector vol- 
ume. To reject dominant downward-going atmospheric 
muon background, only events with at least 300 hits and 
NPE > 3200 are retained and their directions are recon- 
structed with a track hypothesis. For the IC79 sample, a 
log-likelihood (LLH) fit of a muon-track hypothesis based 
on probability distributions of the photon arrival times 
at each DOM is performed . A loose cut on a fit qual- 
ity parameter is applied to filter out events which con- 
tain muons from independent air showers. For the IC86 
sample, a robust regression technique 13, U| is utilized 
to remove hits that have a timing significantly different 
from what is expected from the bulk of the photon signal 
from a muon track. Afterwards, the particle directions 
are reconstructed using the LineFit algorithm Q . These 
algorithms are designed to reconstruct muon track events 
with a zenith angle resolution of 1° or better. For cas- 
cade events which resemble point-like light sources, the 
reconstruction behavior is quite different, and the track 
reconstruction finds nearly arbitrary zenith angles. How- 
ever, there is a tendency toward upward-going and hori- 
zontal directions for the LLH fit and LineFit respectively. 
Since for these directions the NPE cut is lower than for 
downward-going events (see Fig.[5]and Eq.[T|), such events 
can be retained in the final sample even if they would be 
rejected on account of their true direction. 

The final-level selection criteria are based on Monte 
Carlo simulations of background and signal events. They 
are determined by optimizing the NPE-threshold values 
for each of the IC79 and IC86 samples to discover a sig- 
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FIG. 2. Distribution of NPE and reconstructed zenith angle for the IC79 events. The experimental test sample is shown in 
the left panel. The total background (middle panel) and cosmogenic signal neutrino Monte Carlo [f| (right panel) are also 
displayed. The colors show event numbers per livetime of 33.4 days for the IC79 test sample. The signal distribution contains 
all three neutrino flavors. The background distributions include atmospheric muons from the CORSIKA package with the 
SIBYLL high-energy interaction model § assuming pure iron primary composition, conventional atmospheric neutrinos, and 
prompt atmospheric neutrinos. The solid lines in each panel represent the final selection criteria for IC79. 



nal [ID, [l3| from the cosmogenic neutrino model [(| . Fig- 
ure shows the event distributions for the Monte Carlo 
simulations and the 10% experimental IC79 test sample 
corresponding to a livetime of 33.4 days. The solid lines 
in Fig. [5] represent the final selection criteria for IC79 
where events above the lines constitute the final sample. 
The final selection criteria for the IC86 sample are de- 
fined as follows: 



^ 10000 



loginNPE > 




cos 9 < 0.075 
cos $ > 0.075. 

(1) 

The resulting neutrino effective areas, the equivalent area 
at the Earth's surface in which neutrinos are detected 
with 100% efficiency, averaged over the two-year period 
from May 2010 to May 2012 taking into account the dif- 
ferent detector configurations, is shown in Fig. [3l The 
effective areas are given for each neutrino flavor aver- 
aged over 47T solid angle. They show that the analysis 
starts to be sensitive in the energy region around 1 PeV 
with its sensitivity rapidly increasing with energy. The 
effective area is larger for v e than or v T below 10 PeV 
which illustrates the sensitivity of the present analysis to 
cascade events in this energy region. 

The expected numbers of background events in the 
final sample for the 615.9 day livetime from atmo- 
spheric muons and neutrinos from decays of pions 
and kaons are 0.038 ±0.004 (stat) ±8:811 ( s y st ) and 
0.012 ±0.001 (stat) 

±aoo? ( s y s t)i respectively. Compared 
to previous analyses, the utilized atmospheric neutrino 
flux models [l5| accommodate an improved parametriza- 
tion of the primary cosmic-ray spectrum and composi- 
tion which accounts now for the "knee" in the cosmic- 
ray spectrum. Adding prompt atmospheric neutrinos 
from decays of charmed mesons assuming the model 
in 16[ with the improved cosmic-ray spectrum model- 
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FIG. 3. The average neutrino effective area in the present 
analysis for a 4tt isotropic flux, 615.9 days livetime, and the 
IC79 and IC86 string configurations. Exposure of the sample 
used in this analysis is obtained by multiplying the effective 
area with the livetime and 4n solid angle. The sharp-peaked 
structure for i> e at 6.3 PeV is due to the Glashow resonance 

El. 



to 0.082 ±0.004 (stat) 1°;°^ (syst). We note that the 
charmed meson production cross-section used here is un- 
certain to a large degree. Theoretical uncertainties in our 
baseline model [l6[ which uses perturbative-QCD calcu- 
lations are included in the background estimation. In 
addition, potential non-perturbative contributions, such 
as intrinsic charm in nuclei [l7t | or from the gluon den- 
sity at small x, could lead to significantly larger cross 
sections and hence higher prompt neutrino fluxes. The 
IceCube limits on the prompt flux are about a factor of 
four higher than the baseline model [HI]. 

The main systematic uncertainties on the backgrounds 
are from the measurement of NPE and from uncertain- 



5 



"Bert" 



"Ernie" 





FIG. 4. The two observed events from August 2011 (left 
panel) and January 2012 (right panel). Each sphere repre- 
sents a DOM. Colors represent the arrival times of the pho- 
tons where red indicates early and blue late times. The size 
of the spheres is a measure for the recorded number of photo- 
electrons. 



ties in the cosmic-ray flux. Uncertainties in the expected 
number of background events are estimated by varying 
the associated parameters in the simulation. The two 
dominant sources of experimental uncertainties are the 
absolute DOM sensitivity and the optical properties of 
the ice which contribute with (+43%, -26%) and (+0%, 
—42%), respectively. Uncertainties in the cosmic-ray 
flux models are dominated by the primary composition 
(+0%, -37%) and the flux normalization (+19%, -26%). 
The theoretical uncertainty in the neutrino production 
from charm decay [l6| relative to the total background 
is (+13%, —16%). The systematic uncertainties are as- 
sumed to be evenly distributed in the estimated allowed 
range and are summed in quadrature. 

The atmospheric muon and neutrino background 
events are simulated independently. However, at higher 
energies, events induced by downward-going atmospheric 
neutrinos should also contain a significant amount of at- 
mospheric muons produced in the same air shower as 
the neutrino [3]. Since these events are reconstructed 
as downward-going, they are more likely to be rejected 
with the higher NPE cut in this region. Thus, the num- 
ber of simulated atmospheric neutrino background events 
is likely overestimated in the current study. 

After unblinding the 615.9 days of data, we observe two 
events that pass all the selection criteria. The hypothesis 
that the two events are fully explained by atmospheric 
background including the baseline prompt atmospheric 
neutrino flux [lj| has a p- value of 2.9 x 1(T 3 (2.8a). This 
value takes the uncertainties on the expected number of 
background events into account by marginalizing over a 
flat error distribution. Since the prompt component has 
large theoretical uncertainties we have also studied how 
much our baseline prompt component has to be enlarged 
so that the two events can be explained as atmospheric 
neutrinos: obtaining two or more events with a probabil- 
ity of 10% would require a prompt flux that is about 15 
times higher than the central value of our perturbative- 
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FIG. 5. Event distributions for 615.9 days of livetime at fi- 
nal cut level as a function of log 10 NPE. The black points 
represent the experimental data. The error bars on the 
data points show the Feldman- Cousins 68% confidence inter- 
val [2y| . The solid blue line marks the sum of the atmospheric 
muon (dashed blue), conventional atmospheric neutrino (dot- 
ted light green) and the baseline prompt atmospheric neutrino 
(dot-dashed green) background. The error bars on the line 
and the shaded blue region are the statistical and systematic 
uncertainties, respectively. The red line represents the pre- 
diction of a cosmogenic neutrino model (Ahlers et al. [211 ]) 
with the model uncertainty indicated by the shaded region. 
The magenta line represents a power-law flux which follows 
E~ 2 up to an energy of 10 9 GeV with an all-flavor normaliza- 
tion of E 2 4>„ c+Vfl+ „ T = 3.6 x 10~ 8 GeVsr" 1 s" 1 cm" 2 , which 
is the integral upper limit obtained in a previous search in a 
similar energy range [l2| . Signal neutrino model fluxes are 
summed over all neutrino flavors, assuming a flavor ratio of 
v„:v,,:v T — \:\:\. 



QCD model. This contradicts our current limit on the 
prompt flux which would allow for not more than 3.8 
times the central value at 90% C.L. jrsj ]. 

The two events are shown in Fig. S) Both events are 
from the IC86 sample, but would have also passed the se- 
lection criteria of the IC79 sample. The spherical photon 
distributions of the two events are consistent with the 
pattern of Cherenkov photons from particle cascades in- 
duced by neutrino interactions within the IceCube detec- 
tor. There are no indications for photons from in-coming 
or out-going muon or tau tracks. Hence, these events are 
most likely induced by either CC interactions of electron 
neutrinos or NC interactions of electron, muon or tau 
neutrinos. CC interactions of tau neutrinos induce tau 
leptons with mean decay lengths of about 50 m at these 
22j. The primary neutrino interaction and the 



energies 

secondary tau decay initiate separate cascades which in a 
fraction of such events lead to an observable double-peak 
structure in the recorded waveforms. The two events do 
not show a significant indication of such a signature. Fig- 
ure [5] shows the final-cut NPE distributions for the ex- 
perimental data, several signal models and background 
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Monte Carlo simulations. The two events are found near 
the NPE threshold of the analysis. These two events are 
consistent with a previous upper limit by IceCube [12j on 
an unbroken E~ 2 flux (Fig. [5j histogram with downward 
arrows), while a flux corresponding to this upper limit 
predicts about 10 events above the NPE cut. 

For the dedicated reconstruction of the two events, 
a maximum-likelihood fit is utilized. The likelihood is 
the product of the Poisson probabilities to observe the 
recorded number of photo-electrons in a given time in- 
terval and DOM for a cascade hypothesis which depends 
on the interaction vertex, deposited energy and direc- 
tion. Here, the time of the first hit mainly determines 
the vertex position and the recorded NPE plays a dom- 
inant role in estimating the deposited energy. The hit 
information used in the reconstruction is extracted from 
an unfolding procedure of the waveforms. The open cir- 
cles in Fig. [1] indicate the strings closest to the recon- 
structed vertex positions. The reconstructed deposited 
energies of the two cascades are 1.04 PeV and 1.14 PeV, 
respectively, with a combined statistical and systematic 
uncertainty of ±15% each. The errors on the deposited 
energies are obtained by simulating cascade events in the 
vicinity of the reconstructed energies and vertices. The 
study is specifically performed on each event and the 
larger of the two event uncertainties is cited for both 
events. Thus, the error associated with the two events 
is different from errors for other cascade events observed 
in IceCube [23|. Since there is no absolute energy stan- 
dard with adequate precision at these energies, the en- 
ergy scale is derived from simulations based on measured 
ice properties and PMT efficiencies which are assured by 
measurements of atmospheric muons. The main sources 
of the systematic uncertainty on the reconstructed en- 
ergies are the absolute DOM sensitivity and the optical 
properties of the ice The effect of the latter is esti- 
mated to be +9% and —5% and is obtained by varying 
the scattering and absorption coefficients for the photon 
propagation by 10%. The reconstruction algorithm in- 
cludes variations of the scattering and absorption coeffi- 
cients with depth (ice layers) [25;] . The effect of a possible 
azimuthal anisotropy of the ice parameters and a tilt of 
the ice layers on the reconstructed energies is estimated 
to be ±5%. The uncertainty on the absolute DOM sensi- 
tivity is ±10% which scales linearly with the uncertainty 
on the reconstructed deposited energy. The suppression 
of bremsstrahlung and pair production due to the LPM 
effect [26| is negligible in this energy range. The proper- 
ties of the two observed events are summarized in Tab. HI 

The reconstructed energy corresponds to the energy of 
the incoming neutrino if the observed cascade is the result 
of a CC interaction of electron or anti-electron neutrino 
in deep-inelastic scattering, as in this case the total neu- 
trino energy is deposited near the interaction vertex ■ 
On the other hand, NC interactions of neutrinos of any 
flavor or interactions of electron anti-neutrinos via the 



events 



"Bert" 



"Ernie" 



date (GMT) August 8, 2011 January 3, 2012 

NPE 7.0 x 10 4 9.6 x 10 4 

number of recorded DOMs 312 354 
reconstructed deposited 

energy (PeV) 1.04 ±0.16 1.14 ±0.17 

reconstructed z vertex (m) 122 ± 5 25 ± 5 

TABLE I. List of characteristics of the two observed events. 
The vertex z positions are with respect to the center of the 
IceCube detector at a depth of 1948 m. 



Glashow resonance at 6.3 PeV [14| with outgoing lep- 
tons induce cascades which carry only a fraction of the 
neutrino energy. However, given the reconstructed ener- 
gies, the observed cascades are unlikely to originate from 
the Glashow resonance as only about 10% of these in- 
teractions will deposit 1.2 PeV or less in the detector in 
cascade-like signatures. 

The two PeV neutrino events observed in two years 
of data taken with the IceCube neutrino telescope may 
be a first hint of an astrophysical high-energy neutrino 
flux. However, given the yet rather moderate signifi- 
cance of 2.8a with respect to the expected atmospheric 
background and the large uncertainties on its prompt 
component, a firm astrophysical interpretation requires 
more data in combination with analyses in other detec- 
tion channels and energy ranges. 
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